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ABSTRACT 

We present a broad-band spectral analysis of the joint XMM-Newton and Nuclear 
Spectroscopic Telescope Array observational campaign of the narrow-line Seyfert 1 SWIFT 
J2127.4+5654, consisting of 300 ks performed during three XMM-Newton orbits. We detect a 
relativistic broadened iron K a line originating from the innermost regions of the accretion disc 
surrounding the central black hole, from which we infer an intermediate spin of a — 0.58 ' 

The intrinsic spectrum is steep (F = 2.08 ± 0.01) as commonly found in narrow-line Seyfert 
1 galaxies, while the cutoff energy ( E c — 1 08 Iq keV) falls within the range observed in 
broad-line Seyfert 1 galaxies. We measure a low-frequency lag that increases steadily with 
energy, while at high frequencies, there is a clear lag following the shape of the broad Fe K 
emission line. Interestingly, the observed Fe K lag in SWIFT J2 127. 4+5654 is not as broad 
as in other sources that have maximally spinning black holes. The lag amplitude suggests 
a continuum-to-reprocessor distance of about 10-20 r g . These timing results independently 
support an intermediate black hole spin and a compact corona. 

Key words: accretion, accretion discs -galaxies: active -galaxies: individual: SWIFT 
J2127.4+5654- galaxies: Seyfert. 


1 INTRODUCTION 

According to the commonly accepted paradigm, luminous active 
galactic nuclei (AGN) are believed to host a supermassive black 
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hole at their centre, surrounded by a geometrically thin accretion 
disc. The nuclear hard power-law continuum that dominates the 
spectral emission above 2keV is thought to arise in a hot corona 
above the accretion disc, where UV/optical seed photons from the 
disc are Compton scattered towards the X-ray band. This primary 
X-radiation in turn illuminates the disc and it is partly reflected 
towards the observer’s line of sight (two-phase model, Haardt 
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& Maraschi 1991; Haardt, Maraschi & Ghisellini 1994). Such 
a physical process leads to a power-law spectrum extending to 
energies determined by the electron temperature in the hot corona. 
The power-law index is a function of the plasma temperature T and 
optical depth r . 

The X-ray spectra of unobscured AGN is characterized by ubiq- 
uitous and, generally, non-variable features, like the neutral iron 
Kct narrow core and the Compton reflection component (Perola 
et al. 2002; Bianchi et al. 2007). These features can be attributed to 
the reprocessing of the nuclear radiation by distant, neutral matter 
(i.e. the ‘pc-scale torus’ in the framework of standard Unification 
Models; Antonucci 1993). Furthermore, there are several important 
spectral features, like the warm absorber, the soft excess and the 
relativistic component of the iron Kct line, which are present in a 
number of objects (Blustin et al. 2005; Crummy et al. 2006); dif- 
ferent observations of the same object often show these features to 
be variable (Piconcelli et al. 2005; Nandra et al. 2007; de La Calle 
Perez et al. 2010). 

X-ray wavelengths offer the best opportunity in which to inves- 
tigate the physical properties of the primary nuclear source closest 
to the event horizon and to give constraints on its geometry. Recent 
works suggest that most of the emission arises from within ~20r g 
(Reis & Miller 2013) and X-ray microlensing experiments confirm 
that it is compact; in some bright quasars a half-light radius of the 
corona of r s < 6 r g has been measured (Chartas et al. 2009), where 
the gravitational radius is defined as r g = GM/c 2 . Eclipses of the 
X-ray source have also placed constraints on the size of the hard 
X-ray emitting region: r s < 30 r„ in NGC 1365 (Risaliti et al. 2007; 
Maiolino et al. 2010; Brenneman et al. 2013), r s < 1 1 r g in SWIFT 
J2127.4+5654 (Sanfrutos et al. 2013). 

Recently, X-ray reverberation around accreting black holes has 
also revealed that the corona is compact. On short time-scales, 
the variability associated with the primary continuum is found 
to lead the variability associated with the soft excess and the 
broad iron K spectral features. These time delays are short (e.g. 
tens of seconds in 1H0707-465, which has a black hole mass of 
~2 x 10 s Mq; Kara et al. 2013a), suggesting that the light travel 
time between the continuum emitting corona and the reprocessing 
region is small, within ~ 1 0 r g . The first detection of X-ray reverber- 
ation was found between the continuum and the soft excess (Fabian 
et al. 2009). In this work, the soft excess was shown to be associ- 
ated with the broad iron L emission line which peaks at ~0.8 keV. 
Since then, reverberation lags associated with the soft excess and 
the iron K emission line have been discovered (Zoghbi et al. 2012; 
De Marco et al. 20 1 3 ; Kara et al. 20 1 3c), and all reveal that the X-ray 
emission is coming from very small radii. In this paper, we explore 
the frequency-dependent time lags in SWIFT J2127.4+5654 that 
are associated with the soft excess and the broad iron K line. 

In addition to these geometric constraints, cutoff energies in sev- 
eral bright Seyfert galaxies have been measured with hard X-ray 
satellites in the past, such as BeppoSAX (Perola et al. 2002; Dadina 
2007) and INTEGRAL (Panessa et al. 2011; de Rosa et al. 2012; 
Molina et al. 2013). Current measurements of cutoff energies range 
between 50 and 300 keV and require spectra extending above 50 keV 
for better constraints, due to the contribution from complex spectral 
components (such as reflection from the accretion disc and from 
distant material) to the broad-band spectral shape. 

Narrow-line Seyfert 1 galaxies (NLS 1 ) are a class of AGN whose 
spectra show peculiar emission-line and continuum properties. They 
are considered to be Eddington limit accretors (Boroson & Green 
1992) and to systematically host smaller black hole masses than 
broad-line Seyfert 1 galaxies (Komossa & Xu 2007, and references 


therein). The primary defining criteria of these objects are the width 
at half-maximum (FWHM) of the Balmer emission lines in their 
optical spectra (<2000 km s -1 ) and the relative weakness of the 
[Om] emission at L5007 (Osterbrock & Pogge 1985). Past X-ray 
analyses have shown their very variable continua, on time-scales 
of hours (Ponti et al. 2012). The X-ray spectra of these objects 
show very steep slopes, between T = 2. 1 and 2.5 (Boiler, Brandt & 
Fink 1996; Brandt, Fabian & Pounds 1996; Brandt, Mathur & Elvis 
1997). 

The Nuclear Spectroscopic Telescope Array ( NuSTAR ; Harrison 
et al. 2013) is the first high-energy focusing X-ray telescope on 
orbit, ~100 times more sensitive in the 10-80 keV band compared 
to previous observatories covering these energies, enabling the study 
of steep spectrum objects at high energies with high precision. 

SWIFT J2127. 4+5654 (a.k.a. IGR J21277+5656, z = 0.0144) 
was classified as an NLS1 galaxy on the basis of the observed 
FWHM of the Her emission line (~1180kms~ 1 ; Halpem 2006). It 
has been first detected with Swift/BXI in the hard X-rays (Tueller 
et al. 2005). The source was detected in the hard band with 
INTEGRAL! IBIS and an averaged spectrum was obtained by Mal- 
izia et al. (2008) in the 17-100 keV band by summing the available 
on-source exposures (and therefore emission states), resulting in a 
steep spectrum, in the range T = 2.4-3. The authors also discussed 
a single epoch optical spectrum: the NLS 1 classification was con- 
firmed and a black hole mass of 1.5 x 10 7 Mq was inferred. The 
Swift/XKT data were analysed by the same authors together with 
the INTEGRAL! IBIS averaged spectrum and the presence of a spec- 
tral break or of an exponential cutoff in the range of 30-90 keV 
was suggested (Panessa et al. 2011). SWIFT J2127.4+5654 was 
then observed in 2007 with Suzaku for a total net exposure of 
92 ks, a Z.Boi/i'Edd — 0.18 and a T = 2.06 ± 0.03 were measured 
(Miniutti et al. 2009). A relativistically broadened Fe Ka emission 
line was detected, strongly suggesting that SWIFT J2 127.4+5654 
is powered by accretion on to a rotating Kerr black hole, with an 
intermediate spin value of a = 0.6 ± 0.2. This result has been also 
confirmed by Patrick et al. (2011) and in a recent work, using a 
~130 ks long XMM-Newton observation (Sanfrutos et al. 2013). 
The Fe Kcr shape has been also interpreted in terms of reprocessing 
in a Compton-thick disc-wind (Tatum et al. 2012), which would 
require super-Eddington luminosities (Reynolds 2012). 

We present, in the following, results from a simultaneous 
NuSTAR and XMM-Newton observational campaign performed in 
2012 November. Taking advantage of the unique NuSTAR energy 
window, we cover the 0.5-80 keV energy bandwidth. The paper is 
structured as follows: in Section 2, we discuss the joint NuSTAR and 
XMM-Newton observations and data reduction, in Sections 3 and 4 
we present the spectral and lags analyses, respectively. We discuss 
and summarize the physical implications of our results in Sections 5 
and 6. 

2 OBSERVATIONS AND DATA REDUCTION 
2.1 NuSTAR 

NuSTAR (Harrison et al. 2013) observed SWIFT J2127.4+5654 
simultaneously with XMM-Newton with both Focal Plane Module 
A (FPMA) and B (FPMB) starting on 2012 November 4 for a 
total of ~340 ks of elapsed time. The Level 1 data products were 
processed with the NuSTAR Data Analysis Software (nustardas) 
package (v. 1.1.1). Cleaned event files (level 2 data products) were 
produced and calibrated using standard filtering criteria with the 
NUPIPELINE task and the latest calibration files available in the 
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Figure 1. From the top to the bottom, NuSTAR FPMA+B and XMM - 
Newton EPIC-Pn light curves in the 3-80 keV, 2-10 keV and 0.5-2 keV 
energy bands are shown. The ratio between 0.5-2 keV and 2-10 keV EPIC- 
Pn light curves is shown in the bottom panel, straight and dashed lines 
indicate mean and standard deviation, respectively. 

NuSTAR calibration data base (CALDB). Both extraction radii for 
the source and background spectra are 1 .5 arcmin. Custom good 
time intervals files were used to extract the FPMA and FPMB 
spectra simultaneously with the three XMM-Newton orbits. After 


this process, the net exposure times for the three observations were 
77, 74 and 42 ks, for a total signal-to-noise ratio (SNR) of 215, 
235 and 186 for the FPMA, and 209, 234 and 186 for the FPMB, 
between 3 and 80 keV. The background-subtracted count rates are 
0.631 ± 0.003, 0.798 ± 0.004 and 0.889 ± 0.005 counts s“* for 
the FPMA, and 0.604 ± 0.003, 0.785 ± 0.003 and 0.876 ± 0.005 
counts s _1 for the FPMB, between 3 and 80 keV (Fig. 1). 

No significant variability is observed in the ratio between the 
5-10 and 10-60 keV count rates (see Fig. 2) during any of the 
three exposures, suggesting that spectral variability in the NuSTAR 
bandpass during these observations is only weak. 

The three pairs of NuSTAR spectra were binned in order to over- 
sample the instrumental resolution by at least a factor of 2.5 and to 
have an SNR greater than 3<r in each spectral channel. The cross- 
calibration constant between the two instruments is consistent with 
2 per cent in the three orbits. 

2.2 XMM-Newton 

SWIFT J2127.4+5654 was observed by XMM-Newton (Jansen 
et al. 2001) for ~300 ks, starting on 2012 November 4, during 
three consecutive revolutions (OBSID 0693781701, 0693781801 
and 0693781901) with the EPIC CCD cameras, the Pn (Striider 
et al. 2001) and the two MOS (Turner et al. 2001), operated in small 
window and medium filter mode. Data from the MOS detectors 
are not included in our analysis since they suffered strongly from 
photon pileup. The extraction radii and the optimal time cuts for 
flaring particle background were computed with sas 12 (Gabriel 
et al. 2004) via an iterative process which leads to a maximization 
of the SNR, similar to the approach described in Piconcelli et al. 
(2004). The resulting optimal extraction radius is 40 arcsec and the 
background spectra were extracted from source-free circular regions 
with a radius of about 50 arcsec. 

After this process, the net exposure times were 94, 94 and 
50 ks for the EPIC-Pn, with 0.5-10 keV count rates of 
5.411 ± 0.007, 6.498 ± 0.008 and 7.787 ± 0.013 counts s^ 1 for 
the three orbits (Fig. 1), respectively. The observation naturally di- 
vides into low-, medium- and high-flux states. The 2-10 keV flux of 
the source during the third orbit (high-flux state) is consistent with 
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Figure 2. The ratio between 5-10keV and 10-60 keV NuSTAR FPMA light curves (in 3000s bins) is plot versus the time from the start of the NuSTAR 
pointing, straight and dashed lines indicate mean and standard deviation, respectively. 
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Figure 3. Data-to-model ratio of the medium flux state (Orbit 2). EPIC-Pn 
and FPMA spectra are shown in black and red, respectively. The continuum 
is modelled with a T = 1.9 power law (fitted between 2-4 keV and 7.5- 
lOkeV). 


the one observed by Suzaku in 2007 (Miniutti et al. 2009). Spec- 
tra were binned in order to oversample the instrumental resolution 
by at least a factor of 3 and to have no less than 30 counts in each 
background-subtracted spectral channel. This allows the applicabil- 
ity of x 2 statistics. We do not include the 1. 8-2.5 keV energy band 
due to calibration effects associated with the Si/Au edge (Smith, 
Guainazzi & Marinucci 2013). 


3 SPECTRAL ANALYSIS 

The adopted cosmological parameters are H 0 = 70km s 1 Mpc+ 
£2 a = 0.73 and £2 m = 0.27, i.e. the default ones in xspec 12.8.1 
(Amaud 1 996). Errors correspond to the 90 per cent confidence level 
for one interesting parameter (A/ 2 = 2.7), if not stated otherwise. 

We started our analysis fitting the medium-flux state spectrum 
with a power law absorbed with only the Galactic column density 
(7.65 x 10 21 cm~ 2 ; Kalberla et al. 2005) between 2-4 and 7.5- 
lOkeV. In Fig. 3, the broad-band data-to-model ratio is shown (the 
continuum is modelled with a T = 1.89 power law), and both a 
broad iron Key line and a strong Compton hump are clearly visible. 

3.1 EPIC PN spectral analysis 

A phenomenological fit is applied in the 3-10 keV range to the 
three spectra, with a baseline model composed by a power law and 
four Gaussian lines, to reproduce the broad and narrow compo- 
nents of the Fe Key emission line and the Fexxv Key, Fexxvi Key 
lines at 6.7 and 6.966 keV, respectively. We left the normalization 
of the power law free to vary between the three orbits to include 
intrinsic flux variations of the source. The fit is good (x 2 /dof = 
331/321 = 1.03) and no strong residuals are present; results are 
presented in Table 1. Comparing the 2013 observation with the old 
Suzaku one, it is interesting noting that the broad component of the 
iron line (cr = 0 . 22++2 keV) is less prominent and the EW of 
the narrow Fe Kct in the last orbit (the one comparable in flux 
with the Suzaku observation) is consistent with the one reported in 
Miniutti et al. (2009). 

We then simultaneously fitted the 0.5-10 keV EPIC-Pn spec- 
tra from the three orbits with a more physical model, leaving 


Table 1. Best-fitting parameters for the 3-10 keV XMM phe- 
nomenological fit. Asterisks indicate fixed values, energies are in 
keV units, fluxes in 1 0 s ph cm~ 2 s" 1 units and equivalent widths 
for the three orbits are in eV units. 


Id. 

Energy 

Flux 

EW! 

ew 2 

ew 3 

Fe Key [Br ] 

6 37+ 0 - 22 
0.12 

1-3^ 

53+ 55 
J — 25 

43+45 

+5—21 

•jq+40 

•5o_i8 

Fe Key [Nar.] 

r Af'+O.OS 
6 - 46 -0.06 

0.81&J 

31+ 15 
~’ 1 — 15 

27+13 

z/ -13 

y)0 + 12 

Z -5— 12 

Fe xxv Key 

6.7* 

0 5 +0-3 
— 0.3 

20 +1 ° 

zu -10 

1+8 

15+ 

Fe xxvi Key 

6.966* 

0 6 +02 
u -°-0.4 

28^0 

24+ 

21+ 


the normalizations of the primary power law and of the disc 
reflection component free to vary, to take into account the flux 
variations of the source. When we use a model composed of an 
absorbed power law and cold distant reflection the fit is not good 
(X 2 /dof = 849/457 = 1.86) and strong residuals around the iron Kct 
band are present. 

We then added a relativistically blurred reflection component 
arising from an ionized accretion disc and a further absorber at 
the distance of the source (Sanfrutos et al. 2013). We used xillver 
for both cold and ionized reflection (Garcia et al. 2013) and rel- 
conv for relativistic smearing (Dauser et al. 2013). The inclusion 
of a blurred reflection component in the model greatly improves 
the fit (/ 2 /dof = 507/451 = 1.12). The addition of two emission 
lines at 6.7 keV (Fe xxv Key, actually a triplet) and at 6.966 keV 
(Fexxvi Key), with fluxes consistent with the ones reported in Ta- 
ble 1, marginally improves the fit x 2 /dof = 498/449 = 1.11. Re- 
placing the intrinsic neutral absorber with an ionized one, we find 
an upper limit to the ionization state of log £ < —0.54, leading to no 
improvement of the fit. The photon index of the primary continuum 
is steep (T = 2.11 ± 0.01), as commonly found for NLS1 galaxies 
and no strong residuals are present throughout the energy band, 
with the exception of a feature in the spectrum from the third orbit 
around ~0.8keV which is not attributable to any known spectral 
feature (Fig. 4). The excess in the soft energy band, shown in Fig. 3, 
is a combination of the steep primary continuum and the blurred 
ionized reflection, no additional components to the model described 
above are needed. We find an intermediate value for the black hole 
spin (a = 0.5 ± 0.2), in good agreement with the value measured by 


SWIFT J 2127.4+5654: EPIC-Pn best Fit 



1 2 ' ' 5 

Energy (keV) 


Figure 4. EPIC-Pn 0.5-10 keV best fit and residuals with a model where 
blurred relativistic reflection is taken into account, no structured residuals 
are seen. 
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Table 2. Best-fitting parameters for reflection and absorp- 
tion models. Column densities are in 10 21 cm~ 2 units, disc 
inclination angles are in degrees, ionization parameters £ 
are in erg cm s 1 units, Fe abundance is with respect to the 
solar value and cutoff energy E c is in keV units. Power- 
law normalizations are in ph cnr 2 s 1 keV 1 at 1 keV and 
q is the adimensional parameter for the disc emissivity. 
The subscript p indicates that the parameter pegged at its 
minimal allowed value. 



XMM 

XMM+NuSTAR 

/V H 

2.2 ± 0.1 

2.13 ± 0.05 
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00 
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4 
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— 1.7 

6 3 +11 

a 

0.5 ± 0.2 

o d 
+ 1 
oo 
»n 
o 

i 

483} 

49 ± 2 

3 

<20 

<8 

^Fe 

0 54+ 0 06 

0 71+ 0 05 

U - /A — 0.05 

Xefl ( X 10 4 ) 

2.9 ± 0.5 

2.7 ± 0.3 

N 2 

iV refl 

3.7 ± 0.5 

3.5 ± 0.3 

M 3 
yv refl 

3.0 ± 0.5 

3.0 ± 0.4 

X 2 /dof 

498/449 = 1.11 

1733/1566= 1.10 


Suzaku, values for the emissivity and inclination of the disc are also 
consistent with previous analyses (Miniutti et al. 2009). Best-fitting 
parameters are presented in Table 2. When we leave the column 
density of the neutral local absorber and the photon indices of the 
primary continua free to vary between the three spectra, a marginal 
improvement of the fit (x 2 /dof = 481/445 = 1.07) is found, with 
no variations from combined best-fitting parameters. 

3.2 Broad-band spectral analysis 

We then included the six 3-80 keV NuSTAR (FPMA and FPMB ) 
spectra in our analysis, extending the spectral coverage up to 80 keV. 
We left the cross-calibration constants between XMM and NuSTAR 
free to vary. A photon index of T = 2.08 ± 0.01 is found, in 
agreement with previous Suzaku observations (Miniutti et al. 2009). 
We get a best fit x 2 /dof = 1907/1567 = 1.21 and strong, structured 
residuals can be seen above ~20keV (Fig. 5, left-hand panel), 
suggesting the presence of a high energy cutoff. 

Therefore, we substituted the primary power law with the 
cutoffpl model in xspec and tied the values of the cutoff ener- 
gies to the ones of the reflection components. 1 The fit greatly 
improves (Ax 2 = 174, for one additional degree of freedom), 
a x 2 /dof = 1733/1566 = 1.11 is found and no strong residuals 
are present in the whole energy band (Fig. 5, right-hand panel). 
Best-fitting parameters can be found in Table 2 and fluxes, lu- 
minosities and Eddington ratios in Table 3. The cross-calibration 
factors between the Pn and the FPMA, FPMB detectors are 


1 A test version of xillver has been used, where the cutoff energy is a 
variable parameter. 


A'fpma = 1.045 ± 0.007 and .Kfpmb = 1-074 ± 0.007. We mea- 
sure a cutoff energy E c = 1083 jo keV and T = 2.08 ± 0.01: Fig. 6 
(left-hand panel) presents the contour plot of photon index versus 
cutoff energy. In order to estimate the relative strength of the Comp- 
ton hump (arising from both distant and blurred ionized reflection 
components) with respect to the primary continuum, we followed a 
method similar to the one described in Walton et al. (2013) for NGC 
1365, using the relative normalizations of the reflection (modelled 
with pexrav, with the inclination angle and iron abundance fixed to 
best-fitting values) and of the power-law components, above 10 keV. 
We find values of/?! = 1.5 ± 0.2, /?2 = 1.6 ± 0.2 and/? 3 = 1.1 ± 0.2 
for the three sets of spectra, respectively. 

Thanks to the NuSTAR+XMM broad-band view, the degeneracy 
between the contribution from the ionized and the cold reflector can 
be broken, using recent X-ray disc reflection models that include 
relativistic blurring and a direct measure of R. We substituted the 
blurred ionized reflection and primary components in the broad- 
band best-fitting model for relxill (a model that includes both 
relativistic effects and reflection from an accretion disc, Garcia 
et al. 2014). The reflection fraction R has been left to vary between 
the three orbits. The overall fit is good (x 2 /dof= 1756/1570= 1.11) 
and no differences from best-fitting parameters presented in Table 2 
are found. The contribution of the blurred reflection component 
from the disc to the total reflection fraction in the three orbits 
is Rf sc = 1.1 ±0.1, Rf sc = 1.2 ±0.1 and Rf sc = 0.9 ±0.1. The 
contour plot of the cutoff energy E c and the /?? sc parameter is shown 
in Fig. 6, right-hand panel. If we consider a lamp-post geometry and 
use relxilllp (with spin and inclination parameters fixed to their 
best-fitting values) we find a height of the X-ray source of 25 ± 1 0 r g , 
consistent with the estimates discussed in Section 5. 

Once a best fit is obtained we try to reproduce the primary contin- 
uum (modelled with cutoffpl in xspec so far) with a more physical 
model where the electron energy (kT e ) and the coronal optical depth 
(r) can be disentangled, using different geometries. Photon index 
and cutoff energy are fixed to their best-fitting values in the re- 
flection components. We use the comptt model in xspec (Titarchuk 
1994), which assumes that the corona is distributed either in a slab 
or a spherical geometry. In such a model, the soft photon input 
spectrum is a Wien law; we fixed the temperature to 50 eV, ap- 
propriate for Mbh ^ 10 7 Mq. In the case of a slab geometry, the 
fit is good (x 2 /dof =1733/1569 = 1.10) and best-fitting values of 
kT t = 683}} keV and r = 0.353!} 19 are found. When a spherical 
corona is considered, we find a statistically equivalent fit (x 2 /dof = 
1734/1569 = 1 . 10) and best-fitting values are kT e = 533^6 keV and 
r = 1.353 o 67- In both geometries, no significant variations from 
the best-fitting values in Table 2 are found and 3-80 keV fluxes of 
the primary continuum are consistent with the ones presented in 
Table 3. As already pointed out in Brenneman et al. (2014) in the 
case of IC 4329A, the difference in optical depth when two different 
geometries are taken into account is primarily due to the different 
meaning of this parameter in the two geometries: the optical depth 
for a slab geometry is taken vertically, while that for a sphere is 
taken radially (see Titarchuk 1994, for a more detailed description 
of the models). 

4 XMM-NEWTON LAGS 

Finally, in addition to this broad-band spectral analysis, we looked 
for time lags using the XMM-Newton observations. The lags were 
computed using the standard Fourier technique, described in Nowak 
et al. (1999), where the phase lag is found between the Fourier 
transform of two light curves. The phase lag is converted into a 
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Figure 5. Left-hand panel: broad-band fit of the NuSTAR and XMM spectra with a model that does not take into account the high energy cutoff. Clear residuals 
can be seen above —20 keV. Right-hand panel: broad-band XMM +NuSTAR best fit when a cutoff energy is properly modelled. 


Table 3. Best-fitting fluxes and luminosities in different 
energy intervals. Fluxes are in 10 11 erg cm -2 s 1 units 
and are observed. Luminosities are in 10 43 ergs -1 units 
and are corrected for absorption. We calculated Lbol from 
the 2-10 keV luminosity, applying the relation of Marconi 
et al. (2004) for the bolometric corrections (this factor is 
also consistent with the lowest bolometric correction for 
this accretion rate presented in Vasudevan & Fabian 2009). 
A value for the black hole mass of 1 .5 x 10 7 Mg has been 
used (Malizia et al. 2008). 



Orb. 1 

Orb. 2 

Orb. 3 

Fg-lOkeV 

2.40 ± 0.02 

2.87 ± 0.02 

3.32 ±0.02 

ft >-10keV 

1.13 ± 0.01 

1.35 ±0.05 

1.57 ±0.05 

j pow 

^3-80 keV 
F20-100keV 

1.61 ± 0.02 
2.67 ± 0.03 

1.92 ±0.02 
3.21 ±0.03 

2.35 ± 0.02 
3.34 ±0.03 

TgO-lOOkeV 

1.25 ± 0.01 

1.50 ±0.05 

1.57 ±0.05 

^bol / ^Edd 

~0.11 

-0.14 

-0.16 


Confidence contours: Chi-Squared 



90 95 100 105 110 115 120 125 130 

Parameter: HighECut (keV) 


frequency-dependent time lag by dividing by 2nf, where /is the 
temporal frequency. 

We start in usual way by looking for lags associated with the soft 
excess (as in Fabian et al. 2009; De Marco et al. 2013). The left-hand 
panel of Fig. 7 shows the frequency-dependent lag between 0.3-1 
and 1-5 keV. A positive lag indicates a time delay in the hard band, 
which we find at frequencies below 3.5 x 10 -5 Hz. The lag goes to 
zero above this frequency. We then probe the frequency-dependent 
lag between 3-5 and 5-8 keV, where there is less obscuration (right- 
hand panel of Fig. 7). This shows the hard 5-8 keV band lagging at 
higher frequencies, as well (at v < 30 x 10 5 Hz). 

To investigate further, we compute the ‘lag-energy spectrum’, to 
find the average time delays associated with the energy spectrum. 
The zero-point of the lag is dictated by the chosen reference band 
(for this analysis, it was the entire 0.3-1 OkeV, excluding the band 
of interest), and so it is the relative lag between energy bins that 
is meaningful. The left-hand panel of Fig. 8 shows the lag-energy 
spectrum for frequencies less than 3.5 x 10 -5 Hz (where we found 
the 1-5 keV band lagging behind the 0.3-1 keV band). The time 


Confidence contours: Chi-Squared 



Parameter: Ecut (keV) 


Figure 6. Ax 2 = 2.30, 6.17 and 1 1.83 contours (corresponding approximately to 68, 95 and 99.7 per cent confidence levels) for the cutoff energy E c and 
photon index F (left-hand panel) and for the cutoff energy and the blurred ionized reflection fraction in the second orbit (right-hand panel). 
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Frequency (Hz) Frequency (Hz) 

Figure 7. Left-hand panel: the frequency-dependent lag between 0.3-1 and 1-5 keV. A positive lag indicates that the hard band lags behind the soft band. We 
see a hard lag at v < 3.5 x 1 0 s Hz, and zero lag above that frequency. Poisson noise begins to dominate the power spectrum at v ~ 5 x 1 0 4 Hz. (The blue 
hash shows the frequency range explored for the lag-energy spectrum in the left-hand panel of Fig. 8.) Right-hand panel: the frequency-dependent lag between 
the 3-5 and 5-8 keV bands. The hard lag extends up to 3 x 10~ 4 Hz. (The red hash shows the frequency range for the lag-energy spectrum in the right-hand 
panel of Fig. 8.) 




Energy (keV) 


Energy (keV) 


Figure 8. Left-hand panel: the lag-energy spectrum for low frequencies, v < 3.5 x 10~ 5 Hz (corresponding to the frequency range shown in blue hash 
in Fig. 7. The lag increases steadily with energy. Right-hand panel: the lag-energy spectrum for frequencies v = [3.5-30] x HI 5 Hz (the red hash in the 
right-hand panel of Fig. 7. This shows little lag between the continuum and soft excess, and a sharp increase in the lag at >5 keV. 


delay increases steadily with energy, without any obvious spectral 
features, similar to the low-frequency lag-energy spectra in several 
other AGN (Kara et al. 2013b; Walton et al. 2013b). We then probe 
higher frequencies,in the range v = [3.5-30] x 10~ 5 Hz, where the 
5-8 keV band lags the 3-5 keV band (right-hand panel of Fig. 8). 
The high-frequency lag-energy spectrum shows little to no lag be- 


tween the soft excess and the 1-5 keV band, but then jumps up 200 s 
at ~5 keV. The lag peaks at 6-8 keV, and the points from 5-8 keV 
are >4cr above the continuum lag. It does not show the obvious de- 
crease above 8 keV that is often seen in high-frequency lag-energy 
spectra, but statistics are low at this high-frequency, so we cannot 
rule out a downturn of the lag at the highest energies. 
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5 DISCUSSION 

We analysed the NuSTAR+XMM spectrum of the NLS1 galaxy 
SWIFT J2127.4+5654. A relativistic reflection component is found 
in proximity of an intermediate spin Kerr black hole (a = 0.583 q I'), 
confirming past Suzaku results (Miniutti et al. 2009; Patrick et al. 
2011) and providing a very accurate measurement of the intermedi- 
ate black hole spin of this source, thanks to the broad (0.5-80 keV) 
spectral coverage. Fig. 9 shows the different components of the 
model, including the primary continuum, the relativistic blurred re- 
flection of the disc and the emission from distant, neutral matter. 
The contribution from the blurred ionized reflection arising from 
the accretion disc to the total reflection fraction has been measured 
and it is described in Section 3.2. The values obtained are similar 
to those estimated for other type 1 AGN in Walton et al. (2013a). 
Fig. 10 shows the goodness-of-fit as a function of the black hole 
spin: the broader spectral coverage with respect to past analyses 


Current Theoretical Model 



Figure 9. Relativistic reflection model between 0.5 and 1 50 keV. Different 
components such as the relativistic blurred ionized reflection (in orange), 
the primary continuum (in red) and neutral reflection from distant material 
(in blue) can be clearly seen. 



Black hole Spin (a) 


Figure 10. Goodness-of-fit variation as a function of the black hole spin. 
Dashed, horizontal lines indicate Ay 2 values of 1. 4, 9 and 25.7 which 
correspond approximately to 1, 2, 3 and 5cr confidence levels. 


allows us to rule out a maximally spinning black hole with a signif- 
icance >5cr, and a non-rotating Schwarzschild black hole (a = 0) 
is rejected with a significance >3cr. We find that the accretion disc 
parameters (disc emissivity and inclination angle) are in agreement 
with the ones discussed in the past and the steep emissivity (e oc 
+ 63 ) suggests that the Fe emission arises from the innermost re- 
gions of the disc (see Miniutti & Fabian 2004, for a more detailed 
description of the light bending model). Other physical interpreta- 
tions have been recently discussed for the joint XMM and Suzaku 
analysis of Fairall 9 (Lohfink et al. 2012). The intermediate black 
hole spin we measure carries information about the accretion history 
of the source, and a value of 0.6-0. 7 suggests a scenario where the 
black hole mass in SWIFT J2 127.4+5654 is due to a recent black 
hole merger (see Volonteri et al. 2005; Berti & Volonteri 2008, for 
a detailed discussion of the different cases) and Volonteri, Sikora & 
Lasota (2007) for a description of the morphology dependent spin 
evolution. 

The high data quality above lOkeV, thanks to the unique 
NuSTAR spectral coverage, permits us to measure a cutoff en- 
ergy E c = 108+io keV with unprecedented accuracy, marginally 
consistent with prior INTEGRAL results where a cutoff energy of 
49^ keV (Panessa et al. 2011) was measured. Using a Comp- 
tonization model to reproduce the spectral shape of the primary 
continuum, we find values of kT e = 68_ 32 keV and r = 0.35_ 01 g 
with a slab geometry and kT e = 533 2 6 keV and r = 1 .35j/o g 2 with a 
spherical geometry. The power-law continuum is an estimate of the 
power dissipated in the corona, and it is ~57 per cent of the total 
3-80 keV luminosity (L 3 °^ okeV =1-92 x 10 43 ergs _1 for the sec- 
ond orbit, Table 3). This represents ~7.5 per cent of the bolometric 
luminosity of the source (L bo i ~ 2.6 x 10 44 erg s _1 ). 

The electron plasma temperatures agree with the ones found in 
Brenneman et al. (2014) for the broad-line Seyfert 1 galaxy (BLSy 1) 
IC 4329 A, but the lower values for r in SWIFT J2127. 4+5654 sug- 
gest a different geometry in the corona. Indeed, NLS 1 galaxies are 
high rate accreting objects with low black hole masses (Bian & 
Zhao 2004; McLure & Jarvis 2004; Netzer & Trakhtenbrot 2007). 
The lower optical depth might indeed explain why the relativistic 
reflection in this object (~33 per cent of the total 3-80 keV lumi- 
nosity) is more prominent than in IC 4329A, where the power-law 
continuum component is higher, being ~87 per cent of the total 
5-79 keV luminosity. 

SWIFT J2 127.4+5654 shows significant lags over a wide range 
of frequencies. X-ray time lags are largely agreed to be caused by 
two separate processes: by the light travel time between the X-ray 
source and the ionized accretion disc and by mass accretion rate 
fluctuations that get propagated inwards on the viscous time-scale 
and cause the soft coronal emission from large radii to respond be- 
fore the harder coronal emission at smaller radii. When we measure 
the lag at a particular frequency, we are measuring the average time 
lag, and so if both effects are contributing to the lag at a particular 
frequency, it can be difficult to disentangle which process is respon- 
sible for the lag. In sources with a large soft excess that is dominated 
by relativistically blurred emission likes (i.e. 1H0707-495 and IRAS 
13224-3809), we can easily disentangle these two effects by look- 
ing at the lag between ~0.3-l and ~l-4keV. In these sources, 
we find that at low frequencies, the hard band lags the soft (due 
to dominating propagation lags) and at high frequencies, the soft 
band lags the hard (due to greater reverberation between corona and 
accretion disc). For SWIFT J2127.4+5654, the 0.3-1 to 1-5 keV 
frequency-dependent lag shows the typical hard lag of ~ 1000 s at 
frequencies below 3.5 x 10~ 5 Hz (left-hand panel of Fig. 7). It does 
not show an obvious soft reverberation lag, and this likely because 
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the soft excess is not very strong, due to absorption or because of 
low iron abundance. However, we can look for reverberation sig- 
natures at higher energies, where the Fe K line is strong, and there 
is little absorption. The frequency-dependent lag in the right-hand 
panel of Fig. 7 shows the 5-8 keV lagging behind the 3-5 keV band 
up to higher frequencies of ~3 x 10~ 4 Hz. While some of this lag 
is likely due to propagation effects (i.e. below 3.5 x 10~ 5 Hz), 
the high-frequency lag (3.5—30 x 10~ 5 Hz) is likely caused by 
reverberation, where the emission from the Fe K line centroid 
(5-8 keV) lags behind the continuum or red wing of the line 
(3-5 keV). The lag-energy spectra support this hypothesis (Fig. 8). 
The low-frequency lag-energy spectrum below 3.5 x 10~ 5 Hz shows 
a nearly featureless spectrum, increasing with energy (indicative of 
intrinsic propagation lags in the corona), while in the higher fre- 
quencies (v = [3.5-30] x 10~ 5 Hz), we see a different shape. There 
is still a little contribution from the soft excess, but we see a large 
increase in lag at 5 keV, also where we see the red wing of the Fe K 
line. The lag does not decrease above 8 keV, as is often seen in Fe 
K reverberation lags, and this could be an effect of the inclination 
of the disc (Cackett et al. 2014), or just too low statistics to confirm 
a decrease in the lag. 

In principle, the amplitude of the lag reflects the light travel time 
between the primary X-ray source and the accretion disc, and can 
give an estimate of the size of the corona (Wilkins & Fabian 2013). 
However, it is important to take dilution into account before convert- 
ing the time delay to a physical distance. For example, at 6.5 keV, 
the ionized reflection component contributes to ~30 per cent of 
the entire variable emission (i.e. not including the cold reflection 
component), while the continuum contributes to remaining 70 per 
cent. Therefore, the intrinsic lag will be diluted by 70 per cent, 
since 70 per cent of the emission comes from the direct continuum. 
Furthermore, the reference band to which we measure the lag at 
each energy is also affected by dilution. It is composed of roughly 
10 per cent reflection and 90 per cent power law. If the intrinsic lag 
in SWIFT J2127.4+5654 is 1000 s, then, including all the effects 
of dilution, the observed lag at 6.5 keV is ~200 s (which is close to 
what we observe). Given the lower reflection fraction at soft ener- 
gies below 1 keV, the observed lag should be closer to 50 s, which 
is also consistent with what we see in the right-hand panel of Fig. 8. 
Assuming an intrinsic lag of 1000 s puts the X-ray source at a height 
of ~13 r g above the disc, for a black hole mass of 1.5 x 10 7 Mg 
(Malizia et al. 2008). 

The high-frequency lag-energy spectrum (v = [3.5 — 

30] x 10~ 5 Hz) shows that the Fe K lag profile is not very broad. 
Importantly, this timing result is independent of the spectral results, 
which also shows the emission line is not broad. All other sources 
with Fe K reverberation lags have a maximally spinning black hole, 
and show Fe K reverberation lags as low as 3-4 keV. This source 
only shows a reverberation lag at above 5 keV, as expected for a 
larger inner radius. This is very strong evidence for an intermediate 
black hole spin. 

6 SUMMARY 

We presented a NuSTAR+XMM broad-band (0.5-80 keV) spectral 
analysis of the joint observational campaign in 2012 of SWIFT 
J2127.4+5654. 

The main results of this paper can be summarized as follows. 

(i) Thanks to the broader spectral coverage with respect to past 
analyses, a relativistic reflection component is found in proximity of 
an intermediate spin Kerr black hole (a = 0.58^o'|7), a maximally 


spinning black hole can be ruled out with a significance >5cr and 
a non-rotating Schwarzschild black hole is rejected with a signifi- 
cance >3tr; 

(ii) The high data quality above lOkeV allowed us to measure a 
cutoff energy E c = 108^jJ keV with unprecedented precision; 

(iii) The high-frequency lag-energy spectrum shows an Fe K 
reverberation lag. Unlike other maximally spinning black holes that 
have broad Fe K lags, SWIFT J2127.4+5654 has a narrower lag 
profile, which independently suggests an intermediate spin black 
hole. 
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